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SUMMARY 


An  angle-of-arrival  sodar  system  was  designed,  built,  and  tested 
with  the  goal  of  determining  boundary-layer  winds.  The  system  measures 
the  backscattered  signal  induced  in  two  closely  spaced  microphones  on  a 
single  parabolic  receiver  antenna;  the  angle  of  arrival  is  calculated 
from  the  relative  signal  amplitudes.  This  is  the  acoustic  analog  of  the 
amplitude-monopulse  radar  technique.  However,  the  acoustic  system  uses 
distributed  (atmospheric)  targets  and  a  fixed  (not  steerable)  antenna. 

Tests  demonstrated  that  the  system  can  receive  atmospheric  echoes 
and  process  the  analog  signals  to  estimate  angle-of-arrival  (hence, 
layer-averaged  wind)  when  signal-to-noise  ratios  are  adequate.  However, 
the  validity  of  these  wind  estimates  was  not  demonstrated  with  correla¬ 
tive  wind  data.  Digital  processing  techniques  were  implemented  with  the 
goals  of  automatic  wind  calculation,  identification  of  adequate  signal- 
to-noise  ratios,  and  noise  subtraction.  Computer  hardware  limitations 
prevented  achieving  these  goals.  However,  we  believe  that  they  could  be 
achieved  by  using  a  computer  with  larger  memory  capacity. 
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I  INTRODUCTION 


A.  Boundary  Layer  Winds 

A  variety  of  Army  activities  require  information  on  winds  in  the 
planetary  boundary  layer.  For  example,  low-altitude  paradrop  and  bomb¬ 
ing  missions  can  be  critically  affected  by  variable  winds  aloft.  For 
these  applications,  knowledge  of  the  vertically  integrated  boundary- 
layer  wind  (Air  Weather  Service  1971;  Cormier  1975)*  could  greatly 
improve  the  effectiveness  of  operations.  Ballistic  and  rocket  launching 
operations  also  require  information  on  the  vertical  wind  structure. 
Ideally,  data  on  the  wind  effect  over  the  projectile's  total  trajectory 
should  be  available  immediately  before  firing,  so  that  appropriate 
corrections  can  be  applied.  For  current  practical  purposes,  it  would  be 
highly  advantageous  to  provide  wind  data  for  the  layers  in  the  immediate 
vicinity  of  a  weapon  system,  within  one  minute  of  firing,  and  to  heights 
of  about  500  ft  (150  m) .  Such  wind  data  could  be  in  the  form  of  a 
layer-mean  wind  vector.  The  ability  to  make  such  measurements  with  a 
simple  field-transportable  device  could  serve  as  a  point  of  departure 
for  further  developments,  such  as  extension  of  the  technique  to  higher 
levels . 

In  addition  to  the  operational  needs  just  mentioned,  the  Army's 
Automatic  Meteorological  System  (Swingle,  1969)  will  have  a  major  need 
for  small-scale  atmospheric  wind  data.  Such  data  are  required  as  inputs 
to  mesoscale  and  microscale  predictions.  Current  operational,  ground- 
based  Army  equipment  cannot  provide  the  required  winds-aloft  data.  Sig¬ 
nificant  improvements  in  prediction  cannot  be  expected  until  improved 
techniques  become  available  to  measure  the  required  atmospheric  input 
parameters,  including  winds  aloft. 


*References  are  listed  at  the  end  of  the  report. 
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B.  The  Acoustic  Angle-of-Ar rival  Technique 

The  acoustic  angle-of-arrlval  technique  works  by  transmitting  a 
pulse  of  sound  near-vertically  into  the  atmosphere  and  measuring  the  way 
that  winds  transport  the  pulse  in  the  horizontal.  Figure  1  shows  the 
basic  geometry.  For  a  pulse  backscattered  at  height  h,  it  can  be  shown 
(McAllister,  1971a, b;  Georges  and  Clifford,  1972;  Mahoney,  1974;  Mahoney 
et  al.,  1973;  Peters  et  al.,  1978)  that  sound  arriving  back  at  the 
transmitting  antenna  has  its  wavefront  normal  tilted  by  a  angle  9  with 
respect  to  the  vertical,  with  components  given  by 

2 

{9x(h),  6y(h)}  =  —  [u(h) ,  v(h)} 

c  ( 1 ) 

where  9  is  in  radians,  c  is  the  speed  of  sound  (c>>u,v), 

h 

{u(h),  v(h) }  =  —  J  {u(z),  v ( z ) }  dz  , 

h  o  (2) 

and  u(z),  v(z)  are  the  x  and  y  components  of  the  wind  at  height  z. 

Hence,  a  measurement  of  the  angle  of  arrival  {9  (h),  9  ( h ) }  provides  a 

x  y 

measure  of  the  mean  wind  between  the  surface  and  height  h.  The  wind  in 
a  lamina  of  thickness  can  in  principle  be  obtained  by  evaluating  the 
difference  in  angle-of-arrival  measurements;  that  is,  if  we  define 

1  r 

{u(h^i),  v(h„/^0}  =  -  J  { u  (  z )  ,  v  (  z  ) }  dz  , 

&  h  (3) 


we  have  [by  Eqs.  (1)  and  (2)] 


{u(h^\h),  v(h^Ji)}  =  - 


h{6x(h) ,9y(h)} 


-  (h+Z^i)  {9x(h-t^h)  ,  9z(mA)} 


(4) 
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FIGURE  1  BASIC  GEOMETRY  OF  THE  ACOUSTIC  ANGLE-OF-ARRI VAL 
WIND-SENSING  TECHNIQUE 

Propagation  vector  “p  (in  the  ground  frame  of  reference)  is  the  re¬ 
sultant  of  the  wave  front  vector  c  (in  the  moving  frame  of  refer¬ 
ence)  and  the  wind  velocity  vector  Cf  (in  ground  frame  of  refer¬ 
ence).  Propagation  path  shown  assumes  wind  velocity  u  indep¬ 
endent  of  height.  Primed  quantities  refer  to  scattered  waves. 
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The  angle  of  arrival  (0,9)  can  be  determined  by  measuring  either 

x  y 

the  phase  difference  or  the  amplitude  difference  between  signals  induced 

in  closely  spaced  receiver  channels.  The  phase  difference  results  from 

the  fact  that  the  wave  fronts  are  tilted  at  angle  (6,0)  with  respect 

x  y 

to  the  horizontal.  Hence,  for  sound  of  wavelength  x  and  receiver  chan¬ 
nels  spaced  by  distance  (d  ,  d  ),  a  phase  lag  of 

x  y 

2n 

(^x*^v^  (dxsin0x,  dysin©y) 

X  ^  ' 

is  induced  between  adjacent  channels.  The  phase-based  technique  of 
measuring  acoustic  angle  of  arrival  has  been  implemented  by  McAllister 
(1971a, b)  and  by  Peters  et  al.  (1978). 

If  two  microphones  are  mounted  in  the  focal  plane  of  the  same  para¬ 
bolic  dish,  with  each  microphone  offset  from  the  axis  of  the  dish  by  a 
small  distance,  the  receiver  beam  pattern  shown  in  Figure  2  will  result. 
The  sensitivity  functions  R^(9)  and  I^©)  will  each  be  tilted  away  from 
the  axis  by  the  "squint  angle"  qx>  Hence,  signals  arriving  at  angle 
(9^,  0^)  with  respect  to  the  vertical  induce  different  amplitudes  in  the 
two  receiver  channels.  The  amplitude  difference  contains  information  on 

both  the  magnitude  and  the  sign  of  the  angle  of  arrival  (9  ,  9  ).  This 

x  y 

tilted-pattern/amplitude-dif ference  method  of  measuring  angle-of-arrival 
is  the  basis  of  the  monopulse  tracking  technique,  which  has  been  imple¬ 
mented  in  many  microwave  radar  systems  (e.g.,  Rhodes,  1959;  Skolnick, 
1970).  However,  at  the  start  of  this  project  it  had  not,  to  our 
knowledge,  been  implemented  in  an  acoustic  radar  system. 


C.  Relationship  to  the  Acoustic  Doppler  Technique 

Because  the  frequency  of  the  received  acoustic  signal  is  Doppler- 
shifted  by  the  wind  in  the  scattering  volume,  acoustic  radar  wind  meas¬ 
urements  can  also  be  made  by  Doppler  analysis  of  the  received  signal 
frequency.  Many  papers  have  described  the  principles  and  practice  of 
acoustic  Doppler  wind  measurements  (e.g.,  Beran  et  al.,  1974;  Hall  et 
al.,  1975;  Peters  et  al.,  1978),  and  several  commercial  systems  are 


FIGURE  2  ANGLE-OF-ARRIVAL  MEASUREMENT  TERMINOLOGY 
USING  AMPLITUDE  COMPARISON  AND  TILTED 
RECEIVER  CHANNEL  BEAMS 
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available.  At  the  start  of  this  project,  Doppler  acoustic  radar  systems 
described  in  the  literature  used  a  transmit/receive  configuration  with 
three  widely  spaced  antennas  (100  m  apart)  to  measure  both  horizontal 
wind  components  with  good  signal-to-noise  ratio.  An  advantage  of  the 
tristatic  Doppler  system  is  the  ability  to  measure  directly  {u(z),  v(z)} 
profiles  [rather  than  layer  averages  (u(h),  v(h)}]  and  to  obtain  strong 
signals  from  unstable  and  neutral  atmospheric  regions  by  using  scattered 
angles  other  than  near  180°.  (At  scattering  angles  near  180°,  as 
obtained  in  a  monostatic  configuration,  the  strongest  signals  are 
obtained  only  where  the  temperature  profile  is  stable  and  in  the  lower 
parts  of  convective  plumes — see  e.g.,  Hall,  1972;  Wycoff  et  al.,  1973.) 

However,  a  major  weakness  of  the  tristatic  Doppler  system  is  the  need 
for  three  widely  spaced,  mutually  aligned  antennas.  This  need  severely 
limits  the  portability  of  such  systems,  requires  long  setup  times,  and 
restricts  the  number  of  suitable  locations.  In  contrast,  the  acoustic 
angle-of-arrival  technique  can  provide  measurement  of  both  horizontal 
wind  components  (albeit  layer-averaged)  from  a  single  (multichannel) 
antenna.  Thus,  this  angle-of-arrival  project  was  undertaken  with  the 
goal  of  developing  an  acoustic  wind  sensor  that  would  be  more  portable 
and  convenient  than  the  Doppler  systems  described  in  the  literature  up 
to  about  1977. 

During  the  course  of  this  project,  improvements  (made  elsewhere)  in 
Doppler  acoustic  radar  systems  greatly  increased  portability  and  con¬ 
venience.  Therefore,  a  midproject  review  of  the  relative  merits  of 
existing  Doppler  and  potential  angle-of-arrival  (AoA)  systems  was  con¬ 
ducted  in  1979.  The  results  are  summarized  in  Appendix  A.  The  conclu¬ 
sion  was  that  a  system  with  a  single  vertical  AoA  antenna  measuring  pro¬ 
files  of  the  two  horizontal  wind  components,  or  a  single  hybrid  AoA- 
Doppler  antenna  measuring  all  three  components,  would  have  worthwhile 
advantages  over  the  two-  or  three-antenna  monostatic  Doppler  system 
required  to  make  the  same  measurements.  The  primary  advantage  of  the 
AoA  or  hybrid  system  would  be  operational  (smaller  bulk,  better  porta¬ 
bility,  less  expensive  construction  and  maintenance),  although  some 
advantages  in  information  might  also  be  available.  However,  it  was 
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recognized  that  full  achievement  of  this  would  entail  solving  several 
majr -  problems  associated  with  the  AoA  technique.  The  remainder  of  thi 
report  describes  the  SRI  efforts  implement  and  test  the  AoA  technique. 


II  PROCEDURES  AND  RESULTS 


A.  First  Year  (March  1978-February  1979) 

At  the  start  of  this  project,  we  had  planned  to  construct  a 
planar-array  antenna  system  to  make  phase-shift  angle-of-arrival  meas¬ 
urements.  [see  Eq.  (3)].  Such  a  system  would  have  been  similar  to 
those  previously  implemented  by  McAllister  ( 1971a, b)  in  Australia  and 
Peters  et  al .  (1978)  in  Germany. 

However,  before  construction,  we  reviewed  the  amplitude  monopulse 
technique,  its  previous  application  in  microwave  radar  systems,  and  its 
potential  for  use  in  acoustic  wind  sensing.  On  the  basis  of  that 
review,  we  decided  to  implement  the  amplitude  monopulse  technique  using 
a  single  parabolic  dish  antenna  with  the  transmitter  horn  on  axis  and 
receiver  microphones  off  axis.  This  design  contrasts  with  the  speaker- 
array  type  of  antenna  used  by  previous  investigators.  Our  reasons  for 
adopting  the  dish  approach  were: 

•  It  made  maximum  use  of  the  limited  antenna  space  in  the 
trailer-mounted  acoustic  enclosure  used  in  this  project. 

•  Construction  was  simplified  through  the  use  of  fewer  acoustic 
transducers . 

•  Larger  antennas  could  be  assembled  using  no  more  transducers — 
only  a  larger  dish. 

•  Many  microwave  tracking  radars  use  this  design  for  angle-of- 
arrival  measurements  with  good  success. 

Although  the  new  acoustic  design  entailed  some  risk,  we  felt  the  risk 
was  warranted  for  the  reasons  listed.  In  addition,  we  later  realized 
that  the  amplitude-sensing  technique  has  potential  advantages  over  the 
phase-sensing  technique  in  terms  of  noise  subtraction  (see  Appendix  A 
and  the  end  of  this  section). 
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A  dish  antenna  with  two  receiver  channels  (as  shown  in  Figure  3) 
was  constructed,  along  with  transmitter  electronics,  two-channel 
receiver  electronics,  and  a  digital  data  system.  An  outdoor  test  facil¬ 
ity  was  also  constructed;  this  facility  included  a  tower-mounted  movable 
acoustic  source,  a  tower-mounted  u-v-w  anemometer,  and  an  acoustic 
enclosure  for  the  antenna  dish  (on  the  ground).  The  source,  mounted 
above  the  antenna,  could  be  moved  to  +10°  off  axis,  using  a  pulley  con¬ 
trolled  from  the  ground,  to  calibrate  the  system's  response  to  different 
acoustic  angles  of  arrival. 

After  initial  tests,  microphone  preamplifiers  were  added  at  the 
antenna  to  increase  signal-to-noise  ratios.  Although  analog  circuits 
were  assembled  to  measure  the  sum  and  difference  of  signal  amplitudes  in 
the  two  receiver  channels,  these  preliminary  circuits  did  not  perform 
properly,  and  this  approach  was  temporarily  abandoned  in  favor  of  digi¬ 
tal  processing. 

Most  of  the  testing  during  the  first  year  was  accomplished  in  Janu¬ 
ary  and  February  1979  during  a  visit  by  Dr.  Gerhard  Peters  of  the  Max 
Planck  Institute  in  Hamburg.  Dr.  Peters,  who  has  developed  and  tested 
both  Doppler  and  AoA  systems  in  Hamburg  (Peters  et  al.,  1978),  was  a 
paid  consultant  to  this  project. 

The  sensitivity  of  each  receiver  channel  was  increased  by  about  10 
dB  by  modifying  the  microphones  and  placing  them  in  A/A  resonant  tubes. 
Because  the  resulting  system  was  background-limited  by  noise  in  the  SRI 
testing  yard,  no  further  sensitivity  improvements  were  sought  during  the 
test  period. 

The  relative  gain  of  the  two  receiver  channels  was  found  to  be 
unstable  within  about  3  dB.  Although  this  instability  is  unacceptable 
for  long-term  AoA  measurements  using  the  amplitude-comparison  technique, 
it  was  tolerated  for  the  remainder  of  the  first-year  test  period.  Some 
shortcomings  were  found  in  the  narrowband  (commutating)  filter  and  the 
electronic  circuits  for  AoA  detection.  However,  as  both  of  these  func¬ 
tions  were  also  performed  by  software  in  the  data-system  computer, 
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immediate  circuit  improvements  were  not  sought.  Real-time  digital 
methods  for  filtering  and  AoA  measurement  gave  tolerable  results 
throughout  this  test  period;  however,  a  need  for  faster  and  more  effec¬ 
tive  methods  in  the  future  was  recognized. 

The  movable  source  above  the  antenna  was  used  to  measure  the 
receiver  beam  pattern,  R^(0),  for  each  channel  i  (see  Figure  2).  In 
these  first-year  tests,  the  calculated  squint  angles  were 

o 

Ql  =  “92  =  6.7  (6) 


resulting  from  a  microphone  separation  of  6.7  inches  (17  cm)  and  a  dish 
focal  length  of  29  inches  (74  cm).  However,  the  measured  patterns  were 
not  symmetric  about  the  antenna  axis  (8  =  0),  and  neither  pattern  peaked 
at  the  calculated  squint  angle,  +6.7°.  Instead,  the  measured  squint 
angle  was  about  +4°.  The  nonideal  beam  patterns  were  caused  by  reflec¬ 
tions  from  the  test-source  tower  as  well  as  the  transmitter  horn. 

Despite  the  nonideal  beam  patterns,  we  were  able  to  perform  several  use¬ 
ful  tests  with  the  first-year  setup  and  deduce  some  fundamental  diffi¬ 
culties  and  advantages  of  the  acoustic  amplitude  monopulse  technique. 

An  important  difference  between  the  acoustic  monopulse  technique 
for  wind  sensing  and  previous  microwave  radar  raonopulse  tracking  of 
point  targets  is  the  distributed  nature  of  the  target  in  the  acoustic 
wind-sensing  case.  For  the  6-ft  (1.8-ra)  diameter  dish  and  1.0  to  1.6 
kHz  frequency  used  in  our  tests,  the  transmitted  beamwidth  is  initially 
about  13°  (FWHM);  this  width  can  increase  through  turbulent  spreading  as 
the  pulse  propagates  through  the  atmosphere.  Because  of  the  inhomogene¬ 
ous  nature  of  thermal  turbulence  (the  scattering  "target"  in  the  acous¬ 
tic  backscattering  case),  the  source  of  backseat tered  sound  is  not  uni¬ 
formly  distributed  across  the  transmit  beamwidth,  but  may  consist  of  one 
or  more  "bright  spots"  located  anywhere  within  the  beamwidth.  Thus,  the 
acoustic  target,  instead  of  being  located  at  a  unique  angle  0,  is 
characterized  by  a  "brightness  distribution,"  B(8).  The  amplitude  A^  of 
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the  signal-induced  receiver  channel  i  is  the  convolution  of  this  bright¬ 
ness  distribution  with  the  channel  beam  pattern;  hence, 

At  =  /  B(0)  R(e-qi)  d0  ,  (7) 

where  is  the  squint  angle  on  which  the  channel  i  beam  pattern  is  cen¬ 
tered.  Note  that  in  the  case  of  a  point  target  (or  source)  at  0  E(6) 

reduces  to  the  Dirac  delta  function  &(0t),  and  Eq.  (7)  reduces  to 

Aj  =  f  6(et)  R(e-qi)^e  =  R(et-qi)  ,  (8) 

as  expected.  This  is  in  fact  how  functions  RCQ-q^)  were  measured,  using 
the  movable  point  source  (delta  function)  above  the  antenna. 

The  task  of  an  a  .iplitude-based  angle-of-arrival  wind  censor  is  to 
use  the  measured  signal  amplitudes.  A.,  to  infer  something  about  the 
target  brightness  distribution,  B(8).  Fortunately,  it  is  not  necessary 
to  infer  the  complete  function  B(0);  instead  we  try  to  infer  its 
centroid , 


ec  j  B(e)  e  de  / f  B(e)  d«  , 


(9) 


as  being  a  good  measure  of  the  "true"  angle  of  arrival.  Our  tosk  is 
thus  to  construct  some  function  of  the  A.  that  is  unit uelv  related  to 

l 

V 

During  the  January-February  1979  tests  two  such  functions  were 
tried.  The  first, 


A;  -  A2 

f\  -  . 

Aj  +  A2  (10) 
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is  the  usual  difference-over-sum  function  used  in  microwave  monopulse 
tracking  radar  (Rhodes,  1959;  Skolnick,  1970).  In  tracking  radar,  where 
the  antenna  is  steered  using  fj  as  a  feedback  signal,  it  is  necessary 
only  that  fj  ->  0  as  an<*  ^  change  sign  when  ©c  does.  How¬ 

ever,  in  our  application  the  antenna  is  fixed;  hence,  the  exact  shape  of 

f.(0  )  must  be  known,  so  that  f,  can  be  inverted  to  obtain  0  .  We 
1  c  1  c 

therefore  measured  fj  as  a  function  of  ©c ,  the  position  of  the  movable 
acoustic  source. 

The  results  of  these  preliminary  measurements  were  fairly  encourag¬ 
ing,  indicating  that  f^  could  be  inverted  to  yield  ©c  over  a  range  of 
roughly  +4°.  Using  Eq.  (1),  this  corresponds  to  a  wind  speed  range  of 
+11  m/s  or  +25  mi/h.  Unfortunately,  it  was  found  that  f  is  sensitive 
to  the  background  noise  often  present  during  field  operations.  To  see 
this,  note  that,  for  nonnegligible  noise,  Eq.  (7)  must  be  modified  to 
yield 


Ai  *  si  +  Ni  .  (11) 

where  is  the  background  noise  in  channel  i  and  is  the  backscat- 
tered  signal, 

Si  =  f  B(0)  R(0-qi)  d0  .  (12) 

Hence,  with  noise  inputs  present  in  each  receiver  channel, 

Si  -  S2  +  Ni  -  N2 

F1 - 

Si  +  S2  +  Ni  +  N2  (13) 

When  the  tests  at  the  SRI  test  site  were  repeated  for  low  signal-to- 
noise  ratios,  the  effects  of  noise  in  Eq.  (13)  were  evident,  noticeably 
distorting  the  f^  vs  ©c  relationship  that  had  been  obtained  at  high 
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signal-to-noise  ratio.  Thus,  for  nonnegligible  background  noise,  is 
not  a  unique  function  of  ©c» 

The  background-noise  problem  with  f^  could  potentially  be  solved  by 
measuring  in  each  channel  (just  before  pulse  transmission)  and  sub¬ 
tracting  it  from  to  obtain  S^.  However,  fj  has  an  additional  prob¬ 
lem.  Recall  that,  because  of  the  13°  transmit  beamwidth  (FWHM)  and  the 
possible  occurrence  of  "bright  spots"  anywhere  within  the  beam,  ©c  can 
range  from  about  +6°  to  -6°  on  any  given  shot,  even  when  no  wind  is 
blowing.  Because  the  probability  distribution  of  bright  spots  within 
the  beam  should  be  random,  averaging  ©^  over  many  shots  should  yield  the 
mean  wind-induced  ©  .  However,  determination  of  ©  on  a  shot-by-shot 
basis  is  not  possible  using  f^  unless  the  range  of  ©c  over  which  f^  is 
monotonic  exceeds  the  transmit  beamwidth  plus  the  maximum  wind-induced 
angle-of-arrival .  If  this  is  not  the  case,  the  inverse  function  ©c(f^) 
is  not  defined  over  the  possible  range  of  ©  ,  and  the  many-shot  average 


©c  is  not  equal  to  ©c ( f  ^ , 


Because  of  these  problems  with  fj ,  another  function,  f2,  was 


tested.  f2  is  defined  as 


f2  s 


2  2 
A1  a2 

E1  ‘  e2 


(14) 


where  is  the  received  vector  signal  in  channel  i,  i.e., 

E^  =  cos  (wt  -  (ii)  , 

where  inis  the  carrier  frequency  and  4>  is  a  phase  angle.  However,  this 
function  was  abandoned  because  tests  showed  that 

•  Like  f , ,  f2(©c)  was  a^so  nonmonotonic  over  the  range  +6°  (with 
the  6.7  inches  microphone  spacing  that  was  used  in  their  early 
tests),  so  that  the  inverse,  ©c(f2)»  was  undefined. 

•  Noise  can  make  the  instantaneous  denominator  of  f„  zero  or  nega¬ 
tive,  making  very  unstable  at  low  signal-to-noise  ratios. 


A  third  function 


q 

f3  =  - 


a2  -  A1 


2  2A3  -  A2  -  Aj 


(15) 


was  also  defined  during  the  January-February  1979  test  period.  It  was 
not  tested  then  because  f^  requires  a  third  receiving  channel  on  the 
antenna  axis,  f  is  promising  because,  if  each  channel  beam  pattern  can 
be  approximated  by  a  parabola  centered  on  its  squint  angle  q^ ,  f^  is 
equal  to  the  centroid  ©c  of  the  brightness  distribution  B(©),  for  all 
angles  over  which  the  parabolic  approximation  is  valid.  To  show  this, 
we  define  the  parabolic  beam  patterns  as 

2 


Rj  =  C  -  D(9  +  q) 


(16) 


R2  =  C  -  D(©-q) 


(17) 


R3  =  C  -  D© 


(18) 


Substituting  these  into  Eq.  (7)  yields 


A2  -  Ai  =  4Dq9cB 


(19) 


2A3  -  A2  "Aj  *  2Dq  B  , 


(20) 


where  we  have  used  the  definition  Eq.  (9)  and 

B  =  f  B(9)de  . 


(21) 


Substituting  Eqs.  (19)  and  (20)  into  Eq.  (15)  yields 


f3  *  ec  > 


(22) 
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as  desired.  Hence,  is  a  monotonic  (thus  invertible)  function  of  ©c 
for  all  angles  that  (approximately)  satisfy  Eqs.  (16)  through  (18). 
Although  Eqs.  (16)  through  (18)  will  not  be  exact  over  all  angles,  they 
should  be  approximately  correct  over  the  range  where  each  Rj(6)  is 
strong  (i.e.,  the  +6°  beam  width).  Thus,  f^  should  at  least  be  mono¬ 
tonic  over  a  wider  range  of  angles  than  f^  or  and,  after  empirical 
determination,  invertible  for 

All  of  the  January-February  1979  tests  described  so  far  used  the 
movable  elevated  acoustic  source  to  characterize  receiver  response. 
Toward  the  end  of  the  tests,  the  AoA  sodar's  response  to  atmospheric 
echoes  was  tested  and  found  to  be  very  poor.  Thus,  at  the  end  of  the 
first  year,  five  tasks  were  defined  for  the  second  year.  These  were: 

(1)  Improve  performance  on  atmospheric  echoes. 

(2)  Correct  gain  drifts  in  receiver  channels. 

(3)  Improve  digital  filtering  and  amplitude  measurements. 

(4)  Improve  acoustic  enclosure  and  test-source  mounting  setup  to 
eliminate  echoes;  locate  at  quieter  site. 

(5)  Implement  third  receiver  channel  and  test  function  f^. 

B.  Second  Year  (March  1979-February  1980) 

Very  little  technical  effort  was  expended  during  the  first  nine 
months  of  the  second  year  because  funds  were  being  conserved  for  a  major 
effort  that  began  near  the  end  of  the  year  and  continued  into  the  third 
year.  For  purposes  of  continuity,  that  effort  is  described  Section  II- 

C. 


Several  planning  meetings  were  held  during  the  second  year  to 
determine  the  future  direction  of  the  project.  One  of  these  meetings 
was  held  at  the  University  of  Hamburg  while  Dr.  Russell,  the  Principal 
Investigator  on  this  project,  was  visiting  Germany  on  another  project. 
The  results  of  that  meeting  are  described  in  Appendix  A.  As  noted 
there,  it  was  decided  that  the  amplitude-sensing  AoA  technique  had 
enough  potential  advantages  over  acoustic  Doppler  and  phase-sensing  AoA 
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techniques  to  warrant  continuing  the  research,  in  particular  by  attack¬ 
ing  the  five  tasks  listed  at  the  end  of  the  previous  section. 


C.  Third  Year  (March  1980-February  1981) 

Both  electronic  and  acoustic  improvements  were  made  to  optimize 
performance  on  atmospheric  echoes  and  to  correct  the  gain  drift  in  the 
receiver  channels  [Tasks  (1)  and  (2)  of  the  list  in  Section  II-A] .  The 
electronic  improvements  were  made  in  the  preamplifiers,  the  commutating 
filter,  and  the  output  (analog)  channels.  The  preamplifiers  were  com¬ 
pletely  rebuilt  to  reduce  noise,  eliminate  the  signal  gain  drift,  incor¬ 
porate  low-noise  power  supplies,  incorporate  passive  narrowband  filters, 
and  improve  the  input  impedance  matching.  These  improvements,  along 
with  improvements  in  the  microphone  housing,  reduced  the  electronic 
noise  to  well  below  the  background  acoustic  noise,  even  in  quiet 
locations . 


The  commutating  filter  (an  active  narrow-band  filter)  was  rebuilt 
as  a  16-pole  filter.  The  advantages  over  the  previous  four-pole  filter 
are  reduced  sensitivity  to  incoming  signal  phase  shifts  and  simplified 
filtering  of  switching  harmonics  (now  at  16  kHz  for  a  1  kHz  signal). 

The  FWHM  bandwidth  is  30  Hz  at  the  1600  Hz  center  frequency. 


Provisions  were  made  to  record  the  following  output  analog  signals: 


Ramp-gain  amplified  amplitudes,  and  E2 
Aj  and  A2  with  0.1-s  smoothing 

S1  -*2 
*1  +  *2 

-  s2)/(x1  +  a2) 

Trigger  and  timing  data  for  the  digital  system. 


The  block  diagram  for  the  resulting  system  is  shown  in  Figure  3. 
See  Appendix  B  for  a  further  discussion  of  the  analog  and  digital 
circuitry. 
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The  ability  of  this  system  to  detect  atmospheric  echoes  was  demon¬ 
strated  by  comparing  the  output  of  the  SRI  system  with  that  of  a  commer¬ 
cially  available  unit  ( AeroVironment ) .  These  data  are  shown  in  Figure 
A.  The  first,  third,  and  fourth  blocks  of  data  are  from  the  SRI  unit. 
The  second  block  of  data  is  from  the  commercial  unit.  Note  that  the  SRI 
system  readily  detects  echoes  from  both  inversions  and  convective 
plumes.  Although  the  SRI  system  appears  to  show  more  noise  than  the 
AeroVironment  at  high  altitudes,  tests  revealed  that  the  difference  was 
caused  by  an  AeroVironment  ramp  gain  that  increased  less  rapidly  than 
the  desired  linear  function. 

Because  the  SRI  site  used  for  these  tests  was  quite  noisy,  the  sys¬ 
tem  was  moved  to  a  quieter  location  at  Stanford  University.  All  of  the 
following  data  were  taken  at  that  site. 

The  receiver  beam  patterns  were  measured  by  suspending  a  small 
sound  source  21  ft  (6. A  m)  over  the  receiver  dish  and  moving  the  sound 
source  across  the  axis  of  the  receiver  when  no  wind  was  blowing.  This 
sound  source  was  driven  by  the  commutating-filter  oscillator  to  ensure 
that  the  received  signal  was  at  the  peak  of  the  filter  response. 

The  first  measurement  made  with  the  new  system  used  a  single 
transmit/receive  antenna  dish,  with  the  transmitter  horn  mounted  on  axis 
and  the  receiver  microphones  mounted  off  axis  on  either  side  of  the 
horn.  Tests  showed  that  there  were  large  sidelobes  in  the  receiver  pat¬ 
terns,  caused  by  reflections  from  the  walls  of  the  transmitting  horn. 

To  eliminate  these  sidelobes  and  the  associated  distortion  of  f ^ ,  it  was 
necessary  to  separate  the  transmitter  and  receiver.  The  transmitter 
horn  was  mounted  over  a  6-ft  (1.8-m)  dish  and  placed  inside  an  acoustic 
enclosure  beside  the  trailer  containing  the  receiver  dish  and 
microphone . 
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IGURE  4  COMPARATIVE  SODAR  DATA 


The  two  microphones  were  initially  placed  6  inches  (15  cm)  apart  at 
the  correct  focal  distance  from  the  dish.  The  measured  response  pattern 
for  this  configuration  is  shown  in  Figure  5.  Figure  6  is  the  function 

A1  -  a2 
f!  =  - 

Aj  +  A2 

calculated  from  the  data  in  Figure  5.  Note  that  f  is  nearly  linear  for 
1 0 1 1 <4  ,  with  the  inverse  function  given  by  e^ldeg)  =  4.14fj. 

Inspection  of  these  data  indicates  that  the  limits  of  the  linear 
relationship  are  set  by  the  amount  of  overlap  of  the  two  patterns  R^(0) 
and  R  ( 0)  and  the  separation  of  the  minima  between  the  major  lobe  and 
the  first  sidelobe  (see  Figure  5). 

To  test  this  observation,  the  microphone  spacing  was  successively 
set  at  4,  3,  and  2  inches  (10.1,  7.6,  and  5.1  cm)  during  a  calm  period. 
At  each  spacing,  the  sound  source  was  moved  across  the  receiver  beam, 
and  the  function  f^  from  the  analog  output  was  recorded.  The  response 
pattern  for  the  2  inch  (5.1  cm)  spacing  is  shown  in  Figure  7.  Note  that 
the  near-linear,  monotonic  portion  of  f^  is  approximately  14°  wide. 

Thus,  even  if  "bright  spots"  are  present  within  the  +6°  transmit 
beamwidth,  these  bright  spots  should  still  be  within  the  monotonic 
(invertible)  range  of  f^  [cf.  the  discussion  near  Eq.  (6)],  provided  a 
2-inch  (5.1-cm)  spacing  is  used  and  the  wind-induced  angle  of  arrival 
does  not  exceed  about  7°  (i.e.,  wind  speeds  do  not  exceed  21  m/s). 

During  this  phase  of  the  testing,  an  aircraft  flew  over  the  test 
site  and  the  analog  output  of  f  was  recorded.  Although  the  aircraft's 
position  was  not  measured  independently,  the  response  of  the  system 
appeare  to  be  almost  identical  to  that  found  when  the  sound  source  was 
used.  This  showed  that  the  response  was  not  an  artifact  of  the  sound 
source  . 

Figure  8  shows  the  analog  output  data  from  atmospheric  echoes  using 
the  SRI  sodar.  Figure  8(a)  is  the  usual  facsimile  recording  in  which 
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the  magnitude  of  the  return  signal  is  displayed  as  a  function  of  height. 
Note  that  the  data  on  the  left  side  of  this  record  were  taken  with  the 
microphone  spacing  at  6  inches.  The  data  on  the  right  side  were  for  a 
spacing  of  2  inches.  The  intensity  of  the  return  signal  increased 
because  the  microphones  were  closer  to  the  focal  point  and  the  dish  is 
more  efficient  at  that  point. 

Figure  8(b)  shows  the  oscilloscope  record  of  f^  for  nine  successive 
sodar  soundings.  [This  record  was  made  by  putting  a  slowly  increasing 
voltage  on  the  horizontal  axis  so  each  trace  is  displaced  to  the  right 
of  the  previous  one.)  The  fifth  trace  is  a  baseline  trace;  note  that  it 
slopes  to  the  right.  The  soundings  were  made  20  s  apart.  The  calibra¬ 
tion  is  13°/div  or  37  m/s  per  division. 

The  data  in  Figure  8(b)  correspond  to  the  data  at  the  right  edge  of 
Figure  8(a).  Note  that  in  Figure  8(b)  there  is  a  fairly  consistent  peak 
in  fj  at  about  1.9  s,  or  a  height  of  315  m.  Figure  8(a)  also  shows  a 
strong  return  signal  with  a  top  at  315  m.  Hence,  the  peak  in  fj  occurs 
at  a  height  where  signal-to-noise  ratio  is  large  [cf.  the  discussion 
near  Eq.  ( 10) ] . 

The  average  windspeed  parallel  to  the  microphone  (x)  axis  was  cal¬ 
culated  from  trace  7  in  Figure  8(b)  as  follows: 

•  fj  =  1.8  V  (at  1.9  s) 

•  ex  =  2.3°  x  fj  (V) 

•  u  =  1/2  0  (rad)  c 

x 

•  c  =  328  m/s; 

thus , 

9  =4.1° 

x 

and 


u  =  1 1 .9  m/s . 
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Whe-reas  this  is  a  reasonable  value  for  the  mean  wind  x  component  between 
the  ground  and  315  m,  we  were  unable  to  confirm  its  validity,  because 
simultaneous  pilot-balloon  wind  measurements  were  not  made,  and  no 
tower-mounted  anemometers  were  available  at  the  site.  Also,  no  ground- 
based  anemometer  measurements  were  made,  because  the  test  site  was  in 
hilly  terrain  and  it  was  felt  that  ground-level  winds  would  be  poorly 
correlated  with  vertically  averaged  winds. 

Figure  8(b)  also  demonstrates  some  other  difficulties  that  remain 
to  be  overcome  with  this  AoA  sodar.  Note  first  that  f ^ ( 3 1 5  m)  varies 
considerably  from  shot  to  shot,  ranging  from  about  0  in  Trace  8  to  the 
maximum  1.8  V  in  Trace  7.  It  is  doubtful  that  the  vertically  averaged 
wind  (x-component )  between  the  ground  and  315  m  is  changing  from  0  to  12 
m/s  over  the  20-s  shot-to-shot  time  difference.  There  are  two  probable 
reasons  for  the  fluctuations  in  f^(315  m) : 

•  Shot-to-shot  changes  in  the  backscatter  distributions,  B(8) , 
within  the  13°  FWHM  transmit  beam. 

•  Large  shot-to-shot  variations  in  signal-to-noi se  ratio  at  3,' 5  m. 
Both  these  effects  are  discussed  at  greater  length  in  Section 
II-A. 

The  first  effect  could  evidently  be  countered  by  averaging,  for  a 
series  of  shots,  the  8^  value  inferred  from  fj(315  m)  on  each  shot.  The 
resulting  average  can  be  expected  to  have  greatly  reduced  effects  of  the 
varying  brightness  distribut ion,  because  all  measured  fj(315  m)  values 
are  within  the  monotonic  f^  vs  8x  range,  and  the  brightness  distribu¬ 
tion,  B(8),  should  vary  randomly  within  the  beamwidth  from  shot  to  shot. 
However,  shot-to-shot  averaging  will  not  counter  the  second  effect 
because : 

•  Background  noise  itself  can  be  anisotropic  and  correlated  from 
shot  to  shot. 

•  The  fj  vs  0^  relationship  is  a  function  of  noise  whenever  noise 
is  nonnegllgible  compared  to  the  backsca ttered  signal  [see  Eq. 
(12)  and  discussion]. 

Thus,  satisfactory  processing  requires  a  means  of  measuring  background 
noise  in  each  channel  (just  before  pulse  transmission)  and  either 
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subtracting  this  noise  or  excluding  from  angle-of-arrival  calculations 
any  parts  of  the  signal  that  have  poor  signal-to-noise  ratio. 

Small  signal-to-noise  ratios  might  also  explain  some  features  in 
the  vertical  dependence  of  f^(z)  shown  in  Figure  8(b).  Note  the  con¬ 
sistently  small  values  of  f ^ ( 2 9 5  m)  just  below  the  315  m  peak.  If  these 
values  had  been  obtained  at  large  signal-to-noise  ratio  [as  most  of  the 
fj(315m)  values  were],  they  would  imply  very  small  average  wind  between 
the  surface  and  295  m,  and  hence  very  strong  wind  shear  between  295  and 
315  ra — i.e.,  a  jet  at  315  m.  However,  inspection  of  the  sum  output,  Aj 
+  A2,  on  the  facsimile  record  [Figure  8(a)]  reveals  that  the  sum  had 
very  small  signal-to-noise  ratio  at  295  m,  thus,  invalidating  any  infer¬ 
ence  of  0^  from  fj  at  that  height. 

It  is  evident  that  obtaining  useful  wind  estimates  from  this  sodar 
system  requires 

•  Shot-by-shot  noise  identification  (or  subtraction) 

•  Definition  of  signal-to-noise  ratio  as  a  function  of  range 

•  Multi-shot  averaging  of  the  high-signal  inferred  angles-of- 
arrival . 

It  was  for  these  purposes  (plus  digital  calculations  of  A^ ,  fj,  and 

©x)  that  we  included  a  microcomputer-based  digital  daca  system  as  part 
of  the  sodar.  However,  it  turned  out  that  computer  memory  limitations 
pi  evented  implementation  of  the  second  process  listed  above,  and  even 
prevented  satisfactory  determination  of  the  range  at  which  maximum  sig¬ 
nal  occurred  on  each  shot.  Appendix  B  discusses  the  digital  processing 
that  was  implemented  and  the  difficulties  that  were  encountered. 
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Ill  DISCUSSION 


The  goal  of  this  project  has  been  to  produce  and  validate  a  highly 
portable  remote  sensor  of  winds  aloft  using  a  technique  (phase-sensing 
of  acoustic  angle  of  arrival)  that  had  previously  been  described  in  the 
literature  but  not  validated  in  the  practical  sense.  Early  in  the  pro¬ 
ject,  we  modified  the  method  of  approach  to  use  amplitude-sensing  of 
angles  because  of  several  anticipated  advantages  having  mostly  to  do 
with  hardware  (see  Section  II-A). 

Although  the  project  did  not  produce  a  validated,  highly  portable 
sensor,  some  valuable  information  was  gained.  In  particular,  the 
results  show  promise  for  the  amplitude-based  acoustic  angle-of-arrival 
technique,  even  using  the  relatively  simple  difference-over-sum 
function. 

Ai  -  A2 

f!  =  -  . 

At  +  A2 

The  monotonic  range  of  this  function  can  be  made  large  enough  (while 
still  retaining  sufficient  sensitivity  to  angle)  to  span  the  fluctua¬ 
tions  in  angle  of  arrival  that  result  from  the  distributed  and  inhomo¬ 
geneous  nature  of  the  atmospheric  acoustic  target  (i.e.,  the  brightness 
distribution  discussed  in  Sections  II-A  and  II-C).  Thus,  the  three- 
channel  function  f^  described  in  Section  II-A  may  not  be  necessary. 

Moreover,  the  amplitude  AoA  technique  has  an  inherent  advantage 
over  the  phase  AoA  technique  in  terms  of  noise  subtraction.  That  is, 
both  amplitude  AoA  and  phase  AoA  measurements  fail  when  background  noise 
becomes  comparable  to  the  backscattered  signal.  However,  in  the  ampli¬ 
tude  AoA  technique,  noise  effects  can  potentially  be  reduced  by 
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subtracting  the  noise  amplitude  in  each  receiver  channel.  There  is  no 
analogous  phase-noise  subtraction  procedure. 

Demonstrating  the  validity  of  the  amplitude  AoA  technique  requires 
two  advances: 

•  Implementation  of  digital  noise-subtraction,  da ta-select ion ,  and 
data-averaging  techniques. 

•  Successful  comparisons  with  direct  wind  measurements  under  a 
variety  of  conditions. 
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COMPARATIVE  EVALUATION  OF  DOPPLER  AND 
ANGLE-OF-ARRIVAL  SODAR  WIND-SENSING  TECHNIQUES 

In  July  1979,  at  approximately  the  midpoint  of  this  project,  Dr. 
Russell  (Principal  Investigator  for  this  project)  and  Dr.  Gerhard  Peters 
(Principal  Investigator  for  the  University  of  Hamburg's  Doppler  and 
angle-of-ar rival  sodar  projects,  and  a  consultant  to  this  project)  com¬ 
pared  the  perceived  future  for  angle-of-ar ri va  1  (AoA)  measurements  and 
the  future  for  Doppler  measurements.  To  initiate  the  discussion,  they 
prepared  lists  of  AoA  advantages  and  d i sad vantages--bot h  as  perceived  in 
1975  when  SRI's  proposal  to  AR(>  was  written)  and  as  perceived  in  1979. 
Tables  A- 1  and  A-2  show  the  lists.  In  sum,  although  substantial  pro¬ 
gress  with  AoA  measurements  was  made  in  1975-79  (largely  bv  Dr.  Peters' 
group),  much  more  progress  was  made  with  Doppler.  Hence,  the  1975  gap 
between  Doppler  and  AoA  had  widened,  rather  than  narrowed. 

The  review  concluded  by  trying  to  answer  the  following  question: 
"Can  one  expect  to  gain  from  AoA  any  operational  advantage  or  valuable 
information  (on  wind,  turbulence,  or  sodar  technique)  that  cannot  be 
achieved  with  Doppler?"  The  conclusion  was  that  a  svstem  with  a  single 
vertical  AoA  antenna,  measuring  profiles  of  the  two  horizontal  wind  com¬ 
ponents,  would  have  worthwhile  advantages  over  the  two-  or  three-antenna 
monostatir  Doppler  system  required  to  make  the  same  measurements.  The 
primarv  advantage  of  the  AoA  or  hvbrid  svstem  would  be  operational 
(smaller  bulk,  better  portability,  less  expensive  construction  and 
maintenance),  although  some  advantages  in  information  might  also  be 
available  (see  Advantages  3-5  of  Table  A- 1  ).  However,  full  achievement 
of  this  entails  removing  disadvantage  la  of  Table  A-2--that  is,  achiev¬ 
ing  good  AoA  results  in  the  convective  region  below  an  elevated  inver¬ 
sion.  (See  also  advantage  b  of  Exhibit  1  and  disadvantage  1  of  Table 


fluctuations  (Doppler)  as  re¬ 
quired  for  dispersion  studies 


r"v ' 


A-2.)  In  this  connection,  Dr.  Peters  felt  that  the  amplitude  AoA  tech¬ 
nique  pursued  at  SRI  may  have  substantial  advantages  over  the  phase  AoA 
technique  used  up  to  1979  by  the  University  of  Hamburg.  This  is  because 
ambient  acoustic  noise  measured  in  the  various  receiver  channels  is 
always  correlated  to  some  extent,  and  this  produces  some  significant, 
though  spurious,  angle  of  arrival,  especially  when  signal-to-noise  ratio 
is  in  the  neighborhood  of  1.  The  amplitude  AoA  technique  can  remove 
correlated  noise  by  subtraction  in  each  channel,  because  the  amplitude 
measured  by  each  channel  is  a  linear  combination  of  signal  and  noise. 
There  is  no  analogous  "phase  subtraction"  procedure,  because  the  phase 
in  each  channel  is  not  such  a  linear  combination. 

Because  of  these  potential  advantages  of  the  amplitude  AoA  tech¬ 
nique,  the  SRI  research  was  continued  along  the  lines  described  in  Sec¬ 
tion  I I-C  of  the  text. 
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DIGITAL  SIGNAL  PROCESSING 


1 .  Equipment 

The  digital  signal  processing  equipment  is  shown  in  Figure  B-l . 

The  system  is  built  around  a  Digital  Equipment  Corp.  PDP-11/03  computer, 
which  contains  an  LSI-11  microcomputer,  32K  bytes  of  MOS  memory,  and  a 
serial  line  RS-232C  interface.  Additional  capabilities  added  to  the 
system  were  a  floppy  disk  system,  a  nine-track  magnetic  tape  system,  an 
analog-to-digital  converter,  and  a  programmable  real-time  clock. 

Received  signals,  Ej  and  E2 ,  were  sampled  at  the  point  shown  in  Figure  3 
(Section  II-C)  and  input  to  the  ADV11A  analog-to-digital  converter. 
Transmitter  signals  were  used  for  overall  timing  and  sampling  intervals. 
Thus,  these  techniques  are  based  on  the  availability  of  the  transmitter 
signals  for  use  by  the  receiver,  and  therefore,  are  usually  limited  to 
monostatic  radar  configurations. 

2 .  Signal  Modeling 

A  bandpass  signal  can  be  written  as 

E(t)  =  A(t)  cos  [wct  -  ^(t)]  ,  (B-l) 

where  A(t)  is  an  envelope  signal  with  energy  only  in  a  frequency  band 
fw<<fc  (Wozencraft  and  Jacobs,  1963).  <$(t)  is  the  phase  function  which 

specifies  frequency  deviation  from  the  carrier.  Assuming  the  signal 
does  not  have  significant  frequency  modulation  (i.e.,  a  fixed  phase 
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Q-BUS 


FIGURE  B-1  DIGITAL  SIGNAL  PROCESSING  SYSTEM 


delay  i>) ,  trigonometric  identities  allow  E(t)  to  be  decomposed  into  two 
terms : 

E(t)  =  I(t)  cos  (wct)  +  Q(t)  sin  (u^t)  ,  (B-2) 


where 


I ( t )  -  A(t)  cos  i  ,  (B-3) 

Q(t)  2  A(t)  sin  i  . 

I(t)  is  called  the  in-phase  component  and  Q(t)  the  quadrature  component. 
I(t)  and  Q(t)  are  lowpass  signals  with  bandwidth  f  By  synchronously 
sampling  E(t),  the  envelope  of  the  received  signal  can  be  analyzed  with 
filters  at  baseband  [using  I ( t )  and  Q(t)J  instead  of  filters  designed 
near  the  carrier  frequency  (McBride,  1973).  This  greatly  reduces  both 
the  storage  requirement  and  computational  time  needed  to  implement  any 
digital  signal-processing  functions.  The  improvement  is  approximately 
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(fc/fw).  For  the  present  experiments,  f£  a  1600  Hz  and  fy  a  100  Hz 
Thus,  a  reduction  of  16  is  achieved. 


Assume  the  receiver  channel  is  sampled  with  a  period  Tg  =  (l/4)Tc, 

where  T  =  1/f  =  2n/w.  .  The  four  samples  in  time  T  are: 

c  c  c  c 

Es<o>  =  1(0) 

ES(TS)  =  Q(TS) 

ES(2TS)  =  -I(2TS) 

ES(3TS)  =  -Q( 3TS)  .  (B-4) 

This  sampling  process  can  be  represented  in  a  phase  plane  where  one 
revolution  of  the  received  signal  vector  corresponds  to  T^;  E(t)  is  sam¬ 
pled  four  times  during  T^.  Assuming  a  phase  delay  of  zero,  the  four 
sample  points  are  shown  in  Figure  B-2.  A  nonzero  phase  would  rotate  the 
four  points  in  the  plane. 


FIGURE  B-2  PHASE  DIAGRAM  FOR  SAMPLING  OF  RECEIVED  SIGNAL 
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Using  the  series  of  receiver  samples  E  (i.4T  )  and  E  E  (i.4T  +  T  ) 

s  s  s  ss 

(where  i  =  0,1,2,...  counts  the  time  intervals  of  length  4Tg  =  Tc),  the 
in-phase  sequence  I [ i ]  is  defined  as  I(i.4Tg)  =  I(i.Tc),  and  Q [ i ]  = 
Q(i.4Tg  +  T  )  =  Q(i.Tc  +  T  ),  the  quandrature  series.  Because  A(t)  does 

not  appreciably  change  in  T  ,  a  sampled  time  series  of  the  envelope  can 

be  formed  by  the  in-phase  and  quandrature  series  [cf.  Eq .  (B-3)], 

A[i]  =  (I2[  i]  +  Q2 [ i ] ) 1 /2  .  (B-5) 

Thus,  I [ i ]  and  Q[i]  contain  enough  information  to  specify  the  sampled 
envelope  of  the  received  signal. 

As  a  test  of  the  ability  of  synchronous  sampling  to  monitor  only 
signals  about  the  transmitter  carrier  frequency,  a  program  was  written 
which  formed  a  bandpass  filter  centered  at  f  .  The  transfer  function  of 
this  program  is  shown  in  Figure  B-3.  The  response  was  formed  by  averag¬ 
ing  the  sampled  values  from  the  output  of  a  frequency  generator  for  a  1 

s  period.  The  main  response  is  at  the  transmitter  frequency,  f  =  1550 
Hz.  The  response  at  3100  Hz(=2fc)  results  from  aliasing.  Because  the 
receiver  channels  include  a  16-pole  commutating  filter,  this  aliasing 
response  is  removed  before  digital  processing  during  the  experiments. 

By  including  the  in-phase  and  quadrature  measurements  E(2Tg)  and 
E(3Tg),  the  response  at  2fc  can  be  removed  digitally  and  an  overall  nar¬ 
rower  peak  at  can  be  achieved.  The  improved  estimates  are 

I(i  •  4TS)  -  I(i  •  4Ts  +  2TS) 

HD  - - 

2 


Q(i  •  4TS  +  Ts)  -  Q(i  •  4TS  +  3TS) 

Qli]  - - 

2 


(B-6) 


However,  because  these  estimates  require  twice  the  storage  and  because 
the  commutating  filter  minimizes  the  aliasing  effects,  they  were  not 
used  in  the  digital  processing  of  the  receiver  channels. 


FIGURE  B-3  FREQUENCY  RESPONSE  OF  SYNCHRONOUS  DEMODULATED  FILTER 


3 .  Angle-of-Arrival  Estimates 

Table  B-l  describes  the  sequential  steps  of  a  computer  program  that 
calculates  wind  estimates.  An  operator  inputs  the  height  of  the  region 
to  be  studied  along  with  the  width  of  the  transmitter  pulse,  usually  0.1 
s.  The  computer  then  waits  for  a  trigger  from  the  transmitter  indicat¬ 
ing  that  the  pulse  has  been  transmitted.  After  a  time  delay  equivalent 
to  the  propagation  time  of  the  pulse  from  the  transmitter  to  the 
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Table  B-' 


A  REAL-TIME  COMPUTER  PROGRAM  FOR  THE  AMPLITUDE 
TECHNIQUE  OF  ANGLE-OF-ARRIVAL  ESTIMATION 


Input  height  and  transmitter  pulse  duration 

Delay  (2*height)/velocity  after  transmitter  trigger 

Synchronously  sample  receiver  channels  for  duration  of  pulse 

Average  the  in-phase  and  quadrature  components  of  channel 
1  and  2 , 

*1>  *  I2 *  ^2 

2  2  1/2  2  2  1/2 

Calculate  Aj  =  (Ij  +  Qj)  and  =  (I2  +  Q2) 

Ai  -  A2 

Output  function  f  =  - 

1  +  A2 

"Invert"  the  function  into  angle  estimate  using  premeasured 
fj  vs  6  curve 


specified  height  and  back  to  the  receiver,  the  computer  samples  the 
receiver  channels  for  the  duration  of  the  pulse  width.  The  program  then 
computes  the  amplitude  angle-of-ar rival  estimate  using  the  function  f^ 
of  Eq.  (9)  (Barton,  1970;  Carpentier,  1968). 

Amplitude  monopulse  assumes  that  the  two  signals  being  processed 
from  the  target  are  in  phase  but  have  different  amplitudes.  It  is  an 
easier  design  problem  to  obtain  two  amplification  channels  with  the  same 
phase  delay  (especially  when  using  a  narrowband  filter)  than  to  obtain 
two  channels  having  the  same  gain.  However,  the  sequential  sampling  of 
the  two  channels  introduces  a  phase  difference  resulting  from  the  time 
delay  in  the  A/D  converter  (a  single  sample/hold  circuit  with  an  analog 
multiplexer  input),  which  makes  phase  measurements  difficult. 
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The  time  delay  contains  two  parts: 

•  A  fixed  delay  resulting  from  the  execution  of  machine  instruc¬ 
tions  and  the  fixed  conversion  time  of  the  A/D  unit- using  a  suc¬ 
cessive  approximation  technique. 

•  A  variable  delay  resulting  from  the  asynchronous  timing  between 
the  computer  clock  and  the  A/D  clock. 

The  fixed  time  delay  is  54  ys,  an  equivalent  phase  delay  of  31°  at  the 
transmitter  frequency.  The  variable  time  delay  ranges  from  (  to  4  ys,  a 
maximum  phase  delay  of  2.3°. 

The  fixed  and  variable  delay  is  of  little  consequence  with  ampli¬ 
tude  measurements,  because  it  is  assumed  that  the  amplitude  does  not 
change  appreciably  during  this  short  time.  Because  part  of  the  delay  is 
fixed,  it  is  possible  for  a  phase  measurement  program  to  compensate  for 
this  delay.  However,  with  a  phase  measurement  system,  the  error  result¬ 
ing  from  the  variable  phase  delay  could  be  significant.  Thus,  it  is 
important  for  phase  measurement  systems  to  use  A/D  converters  with  mul¬ 
tichannel  sample/hold  circuits,  which  simultaneously  make  analog  meas¬ 
urements  and  then  are  sequentially  converted  into  digital  values. 

During  the  execution  of  the  amplitude  measurement  program,  an 
operator  had  to  input  the  height  of  the  desired  region.  This  was  usu¬ 
ally  the  height  of  an  inversion  layer  which  was  being  displayed  by  a 
nearby  facsimile  recorder  using  the  analog  sum  of  both  receiver  channel 
amplitudes.  Difficulty  in  reading  the  correct  height  of  inversions  led 
to  a  preface  segment  in  the  program  which  monitored  the  strengths  of 
echoes  for  a  single  transmitter  pulse.  The  height  of  the  region  of  max¬ 
imum  return  was  output  to  the  operator  and  also  used  as  input  by  the 
angle-of-arrival  program  during  the  next  transmitter  pulse.  Memory  lim¬ 
itations  allowed  only  sampling  every  50  meters  in  the  vertical,  thus  it 
was  possible  to  miss  the  area  of  maximum  return.  This  caused  a  large 
variation  in  the  angle-of-arrival  estimates  due  to  weak  signal  returns. 
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I)r.  Edward  M.  Liston,  Senior  Engineer 
Dr.  Stephen  A.  DeLateur,  Research  Engineer 
Dr.  Roger  S.  Vickers,  Senior  Physicist 
Dr.  Gerhard  Peters  (consultant) 

Mr.  Henry  D.  Olson,  Senior  Research  Engineer 

Dr.  Norman  I.F.  Chang,  Senior  Research  Engineer 

Mr.  Walter  Jimison,  Research  Engineer 

Mr.  William  Dyer,  Engineering  Assistant 

Mr.  William  Evans,  Staff  Scientist 

Mr.  Ronald  Presnell,  Senior  Research  Engineer 

Dr.  Vi.ent  Salmon,  Staff  Scientist 
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this  pro  j er t . 
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